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Abstract. We start a systematic study of the Lemaitre- 
Tolman-Bondi (LTB) model as applied to the large scale 
structure and its evolution. Here we study three possible 
initial conditions of the LTB models which are asymp- 
totically FRW at large scales: bang time, fractal density 
(with fractal dimension D = 2), and velocity law. Any 
two of these determine the third one. Fractal density and 
simultaneous bang time provide a quantitative estimate 
for the scale beyond which the deflection from the linear 
Hubble law is small. This border may be identified with 
the zero-velocity surface. For fractal density and linear 
Hubble law it is shown that the bang time is necessarily 
non-simultaneous. 
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1. Introduction 

Lemaitre-Tolman-Bondi (LTB) models are exact solu- 
tions of Einstein's equations for 1) spherical symmetry, 
2) pressureless matter (dust) and 3) motion with no par- 
ticle layers inte rsectin g. Originall y stu died by Lemaitre 
(|1933|) , Tolman (|l934| ) and Bondi ( jl947|) , these models are 
the simplest generalization of the Fricdmann-Robertson- 
Walker (FRW) models with a non-zero density gradient. 
At least three cosmological applications of LTB mod- 
els have been discussed. The first relates to the evolution 
of inhomogeneities and peculiar velocities in an expanding 
universe. Important results have been derived under the 
assumption of "unique bang time" (Olson & Silk 1979), 
when every mass shell has been simultaneusly created. So, 
Silk & Wilson p79a| , pJ79b| ), Olson & Silk @jJ7§ ), and 



Olson & Stricland ( 1990 ) studied the formation of galactic 
clusters from small density and velocity perturbations (im- 
plicitly, non-simultaneous bang is used in Silk & Wilson 
(1979b)). It was shown that after a long time the initial 
conditions are forgotten and a universal density profile is 
formed. The LTB model has been applied to the determi- 
nation of the mass density parameter f2 from the local 
peculiar velocity field (Silk 1974), and of the mass of the 
Virgo cluster (H offma n ct al. 1980; Tully & Shaya 1984; 
Teerikorpi et al. [1992; Ekholm et al. 1999). 



The second area of cosmological a pplica tion was devel- 
oped b y Elli s et al. (1985), Maartens (|l996|) and Mustapha 
et al. ( 1998 ). They use the FRW model for small (observ- 
able) scales and assume the LTB model for large scales. 

The third area is the modeling of fractal matter dis- 
tributions within general relativity (Bo nnor [1972 ; Ribeiro 
19924 |l992b| , [1993[ Humphreys et al. [l998b| ; Matravers, 
1998). This has gained impact from redshift surveys re- 
vealing fractality with the fractal dimension D « 2 in the 
space distribution of galaxies up to distances of 100/i -1 
Mpc (h = i7o/100fcrns _1 Afpc _1 ) (see Baryshev et al. 
1994; Sylos Labini et al. 1998), and confirming the scale 
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invariant de Vaucouleurs (|1970| ) law. This leads to a new 
application of the LTB models, pointed out by Bonnor 
(1972), where the fractal structure is treated as spherically 
symmetrical inhomogeneities centered in every galaxy. 

Baryshev et al. ( |1998[ ) showed that the linear pertur- 
bation theory for the growth of density fluctuations leads 
to a non-linear Hubble law if all matter is in the fractals. 
Then the observed linear Hubble law (at scales < lOO/i -1 
Mpc) requires a very low background (FRW) density, 
il < 10~ 3 . The exact LTB calculations by Humphreys 
et al. ( |1998b[ ) and Matravers ( 1998J ) gave the same con- 
clusions. Another explanation for the linear Hubble law 
within the fractal structure (the "Hubble-de Vaucouleurs 
paradox") was proposed by Baryshev et al. (1998): uni- 
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form dark matter with a very high density. Fl The present 
paper shows a third way to make the linear Hubble law, 
by abandoning the assumption of a unique bang time. 

There are two ways to parameterize LTB models. The 
first, introduced by Tolman and Bondi, is called the 3 + 1 
approach . The second uses observational coordinates (El- 
lis et al. |1985| ). As underlined by Matravers (1998), the 
3 + 1 coordinate approach, adopted in the present paper, 
provides a physical interpretation of the evolution of the 
universe in co-moving coordinates. 

Our series of papers is concerned with the 1st and 3rd 
of the above areas and their theoretical basis. In this paper 
we study initial conditions and compatibility of density 
and velocity laws. In Paper II we add the A-term to the 
LTB model, and in Paper III address conceptual problems. 

In Sect. 0we review the LTB models in co- moving co- 
ordinates. In Sect. |we study conditions for the existence 
of the zero-velocity surface in the closed central part of 
the LTB model. We find two criteria to check out if the 
zero-velocity surface exists. In Sect. we show the ve- 
locity deflection from the Hubble law within LTB models 
with unique bang time and density law with fractal dimen- 
sion D ss 2, apply the criterium for the existence of the 
zero-velocity surface, and show how its location depends 
on initial conditions and cosmological density Slo- We also 
calculate the bang time function which reproduces the lin- 
ear Hubble law at all scales within the inhomogeneities 
described by LTB models at the present cosmic epoch. 

2. The LTB models in co-moving coordinates 

The LTB models are the simplest exact, nonlinear, in- 
homogeneous, nonstationary, isotropic dust models in 
general relativity. In co-moving (Lagrangian) and syn- 
chronous coordinates r, t the metric is: 



ds 2 (r, t) = c 2 dt 2 - e A(r '* } dr 2 - R 2 (r, t)dtt 2 



(1) 



and the energy-momentum tensor has the form: T^ = 
diag(c 2 p(r, t), 0, 0, 0). Here c is the speed of light, dVL 2 = 
d9 2 + sin 8 d<fi 2 , and R(r,t) is an Euler coordinate. Such 
metrics and energy-momentum tensor, together with Ein- 
stein's equations, produce a set of inhomogeneous cosmo- 
logical models, generally with time- and space-dependent 
curvature, which are described by the following system 
(Tolman 1934): 
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Sandage et al. (1972) were the first to note the strange co- 



existence of the linear Hubb le law and local inhomogeneities. 
Recently, Teerikorpi et al. (1998) showed from Tully-Fisher 
distances that the all-sky average number density decreases 
as predicted by a fractal with dimension ~ 2.2, from 1 to 
100ft. - 1 Mpc, where (and beyond which) the Hubble law is seen. 
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T^ = ^—[2R'-XR'j=0, 
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where ' = -§- and ' = „? ,, 

or d(ct) 



defines the function A(r, t): 
e A ^=i?' 2 (r,t)// 2 (r). 



= 0, (4) 

(5) 

(6) 
Einstein's equation Eq.(J5|) 

(7) 



f(r) is one of the undetermined functions of the models 
Here ' = # and ' = d 



777777. Via Eq.ffl), the T^-component 
: motion 

(8) 



reduces to the equation of motion 

2 R(r, t)R(r, t) + R 2 {r, t) + 1 - f 2 (r) = 0. 

The first and second integrals of Eq. (pi) are 



R 2 (r,t)= f 2 (r)-l + 



R(r,t) 

±t + t R {r)= f 



2R(r,t)' 



dR 



R(r,t ) 



c 2 (f 2 (r)-l) + ^' 



0) 



(10) 



where F(r) and tfj(r) are the second and third undeter- 
mined functions of the models. tn(r)tn(r) is the bang time 
iii(r). Equation (|2|) gives for the density: 



^ISM-^ 



dr 2 R 2 (r,t)^± 



(11) 



Eq. ([LI]) shows that the models become singular by two dif- 
ferent causes, defined by R(r, r) = (which corresponds 
to the bang time tfi(r)) and R'{r,r) — (which corre- 
sponds to the layer intersection time i^'(r)). The bang 
time (Silk & Wilson 1979a) and layer intersection time 
function (for flat LTB model see Gromov 1997), respec- 
tively, correspond to these singularities. The gravitating 
mass M grav and the invariant mass M; nv of the dust are 
defined by the energy-momentum tensor: 

R(r,t) 

M grav (R(r, t)) = -£ f T °(y) y 2 dy = M grav {r) (12) 



and 



M inv (r) 



47T 



T§(x,t) \f-g{x, t)dx, 



(13) 
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R' R 2 



where y/—g - — - 
of T? and Eq.© 



M grav (r) 



M inv (r) 



Substituting p from the definition 
into Eq.© and Eq.(|l|) we obtain 
(c 2 /4G)F(r), andF(0)=0. So, 



dM, 



grav(x) 



m 



M' 
f( r ) = _21^ 

n> ML' 



(14) 



One of the interpretations of the function f(r), given by 
Bondi (1947), relates it to the curvature and components 
of the Einstein tensor: 
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The space curvature is equal to zero if and only if / = 1. 
The LTB models are defined up to some transforma- 
tion of the co-moving coordinate if> : r — > f , which de- 
creases the number of undetermined functions from 3 to 2 



( Just 1196$ Just & Kraus |l962j ). See also Hellaby & Lake 
( |1985| ) and Hellaby (1987) for interesting examples. These 
two functions should be chosen from the set: 



tR, tRI, PO, /, Rq, Rq- 



(15) 



The transformation is not unique and may be chosen 
according to the character of the problem (see Gromov 

The transformation ifi is time independent, so it can be 
used to fix one of the functions from the set Eq.([15J). In 
this case one speaks about a parametrization of the LTB 
model. We give two examples. The first is often used: 



R{r, 0) = r. 



(16) 



This implies that F(r) is defined by the initial den- 
sity profile p(r, 0). It was fully studied by Liu ( |l990a| , 
1990b|Jl99l| ), and used e. g. by Ribeiro (1992a, 1992b, 1993) 
and Gongalves & Moss ( 1997 ). An alternative transforma- 
tion is: 



M inv (r)G/c 2 =r. 



(17) 



I t was first u sed by Eardley ( 1974 ) and studied by Gromov 
( |1996| , |1997| ). 

The bang time ta is used as one of the initial conditions 
by Silk & Wilson (1979a), Olson & Silk (1979), Olson & 
Stricland (1990). 

3. Zero-velocity surface in the LTB model 

We use two different sets of initial conditions: 

A) bang time tn(r) and initial density profile po(r), and 

B) initial density po(r) and velocity Ro(r) profile. 

They are interconnected as may be seen from a simple 
analogy of an apple dropping from an apple tree. If the 
initial position and velocity, and an equation of motion are 



given, we can calculate the time it will reach the ground. 
For the LTB models it is the same. If we start with tR(r) 
and po(R) as initial conditions, we can calculate the ve- 
locity profile R(r) for t = 0: since a particle must reach 
the centre by the time tR{r), it must have a predefined ve- 
locity at t = 0. Since the density profile is also given, the 
gravitational potential becomes fixed by the same initial 
conditions. As from each two functions the third may be 
derived, all three cannot be selected freely. 

3.1. Dimensionless equations 

Now we restate the models in terms of dimensionless quan- 
tities. We use the following characteristic values: 



l<3 — cto, to — 



H a 



fio = 



p(oo) 



8 7rG Pcr = 3i/ 2 M =~p(oo)Jg, 



and dimensionless variables: 
,. R r 



h 



m 



p(R) 

p(oo)' 



Kv) 



M„ 



Mo 



8{x) x 2 dx, 



where Iq and to are the characteristic length and time, 
il is the density parameter of the FRW background, H 
is the value of the Hubble parameter at the moment of 
the initial conditions, po{r) is the dust density, p cr is the 
critical density, Mq is the characteristic mass, £ and r/ are 
dimensionless Euler and Lagrangian coordinates, S(rj) is 
the dimensionless density, and p(r)) is the dimensionless 
mass of the dust. In terms of these quantities, the bang 
time is written as t^(tj). The index £ reminds us that the 
bang time is the time required for the particle to come 
from its initial position £o to £ = 0. 

In this section we show how one can use the first and 
second integrals, Eq. (H) and Eq. ( |l0| ) , to calculate the zero- 
velocity surface for LTB models if we relax the often used 
assumption of unique bang time. 

We use an effective dimensionless gravitating mass p* 



p* = Q p = 2M grav G/c 2 l . 



(18) 



In terms of dimensionless variables the first integral of the 
equation of motion Eq. M) becomes 



i 2 (v,r) = f( V )-l + 



/'" 



tfary 



(19) 



The form of the second integral Eq.(|lO|) depends on the 
sign of f 2 (rj) — 1. For f 2 (r]) — 1 < (closed models): 



M*W 



-1/2 



A 3 / 2 



( arcsiny A 



±T + T i {r]) = 

1 



AVT 



(20) 
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for f 2 (rj) - 1 = (flat models): 



i 



3/2 _ ^3/2 



C±o^; 



(21) 



and for f 2 (r]) — 1 > (open models): 



±T + r i (i 1 ) = 



A** fa) 



-1/2 



(-A) 3 / 2 
In these formulas: 

c(i) 



-AVl~4 



(22) 



5. ,2. TTie closed and open models with arbitrary bang time 

In this and the remaining sections we concentrate on ini- 
tial conditions. The co-moving coordinate is chosen as £o- 
By solving Eq.(fl9|) for 1 — f 2 and substituting this into 
Eq. ( po|) , the expression for bang time r^ of the closed and 
open models may be rewritten in the form (taking into 
account that for a fixed moment, £ may be used as a co- 
moving coordinate): 



n 




B = _-° > o, 



and for closed models < B < 1 

*(B) = * cl (B) = "■rcsiny/T=B-y/T=By/B 

while for open models B > 1 

*(B) = *°v(B) = -*rcsinh^B—l + jBVB=l 

1 J ~ [ > (B-l) 3 / 2 

Eq.(p4|) allows us to rewrite equation Eq. (|l9|) as 

B = (/ 2 -l)% + l. 

I 1 



(24) 



(25) 



(26) 



(27) 



It follows from Eq.(B4h that B = corresponds to the 
following set of initial conditions: if £o = 0, gradS(^Q — 
0) = 0, when 



lim£~ lim ( Co/Co ) > 0; 

so— h) Co— »0 



(28) 



whereas in the case of £o 7^ 0, B — implies £o = 0. The 
velocity of the particle is equal to zero at the boundary 
£zv, see Figs. 1-3. For both cases f 2 = 1 — /i*/£o > for 
_B = 0, which implies the inequality £o > M* f° r B = 0. 
Note that this restricts the kind of particular TB model in 
which the nonequality may be satisfied: because f 2 > 0, 
it follows that f 2 — 1 < 0. So, B = may be satisfied only 
in the closed model. 

The limit B — > 1 corresponds to / — > 1, so that open 
and closed models both have a common limit which coin- 
cides with the flat model: ^ fl (B) = §. 



Olson & Silk (1979) defined the boundary between 
open and closed TB models with a unique bang time as a 
place where / — 1. We also postulate a set of initial con- 
ditions (e.g. fractal density and Hubble law), which pro- 
duce the following sequence of particular models: a closed 
domain of the model which has a position around a cen- 
tre ("core") and open model farther out from the centre 
(" shell" ) . The two domains are separated by the flat one 
located on the surface where f 2 = 1. For the closed "core" 
< B < 1, so, from Eq.(f|) and Eq.@ it follows that 



(23) < m cl (B) < arcsin(l) « 1.57, 



(29) 



where | corresponds to the boundary of the closed model, 
the flat model, (this boundary we denote by £/;) and arc- 
sin(l) corresponds to the zero- velocity surface £,zv- 

For arbitrary bang time, the zero-velocity surface is 
defined by the solution of Eq. (pG) : 



!£3 

T£ = \l — arcsin(l). 



(30) 



The solution of Eq.(|3(j) may be real or complex depend- 
ing on initial conditions, i.e. bang time and density profile. 
A complex solution means that the zero-velocity surface 
lies in the centre of symmetry. If the solution is real (and 
positive) then the zero-velocity surface is found at a fi- 
nite distance from the centre and separates the collapsing 
region of the closed part of the model from the expand- 
ing region. In the domain < £o < £zv we can introduce 
some other LTB model corresponding to the initial density 



profile, also closed. Humphreys et al. (1998a) first demon- 
strated how to construct the LTB model for that domain. 
The above approach utilizes the coordinate criterium 
for the existence of a central collapsing domain. Using 
Eq.(E4h, we can also define a second form for this cri- 
terium, the mass criterium. From Eq.(J2J) it follows that 
the two limits of function *B(B) correspond to two char- 
acteristic masses. The mass jjlzv ', 



Hzvfo) = (l-57/r € ) 2 Zl 



(31) 



corresponds to the low limit of radial Euler coordinate 
£,zv ■ This designates the beginning of the domain £o > 
£zv from which all particles can collapse at time T£ . Sim- 
ilarly, the characteristic mass corresponding to the flat 
model (or to the upper boundary of the closed model, 
which is the same thing) has the form 



M&) = (2/3t c ) 2 & 



(32) 



This criterium can be stated as follows: if the graph of 
the mass, corresponding to a given initial density profile, 
intersects the graph of /izy, then £zv > 0- 
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3.3. The flat LTB model 

We now turn to the simplest initial condition, / = 1. If 
T£ = const, the flat LTB model reduces to the flat FRW 
model. As shown by Gromov (1997), for the flat LTB 
model the bang time may be given in the form: 



t «(m) 



\ 



M 



dy 

po{vY 



(33) 



which immediately implies that po(£o) = const for simul- 
taneous bang time. In any other case, po(Co) 7^ const and 
the bang time is not constant. As was shown above, £zv 
may not be equal to zero (and the bang is not simultane- 
ous) if and only if the LTB model is closed, so the domain 
of definition of the flat model is the whole region £o > 0. 

4. LTB models for a fractal density distribution 

In this section we study the LTB models with initial con- 
ditions given by the fractal density profile and Hubble law. 

4-.1. Fractal density distribution and the Hubble law: 
Hubble-de Vaucouleurs paradox 

Two fundamental empirical laws have been established 
from extragalactic data. First, there is the power law 
density-distance relation (de Vaucouleurs law) which cor- 
responds to fractal struture (Mandelbrot 1982 ) with frac- 
tal dimension D rj 2 up to the depth of available cat- 
alogs sw 100 h- 1 Mpc. (Sylos Labini et al. 1998). Sec- 
ond, Cepheids, TF-distance indicator and Type la super- 
novae confirm the linearity of Hubble's redshift-distance 
law within the same distances where the fractality exists. 

Baryshev et al. (1998) emphasized that the linear 
redshift-distance relation inside the fractal (inhomoge- 
neous) matter distribution creates the so-called Hubble-de 
Vaucouleurs (HdeV) paradox. It means that the interpre- 
tation of the Hubble law within FRW cosmological mod- 
els as a consequence of a homogeneous galaxy distribution 
disagrees with new data on the galaxy distribution for a 
scale interval from 1 to 100 Mpc. We emphasize that the 
paradox exist for small distance scales (up to 100 Mpc), 
i.e. for redshift s less then 0.03. Hence the arguments of 
Abdalla et al ( 1999 ) on essential relativistic corrections 
do not explain the paradox. 

Two solutions of the HdeV paradox are previously 
known. The first one (Baryshev et al. 1998) is based on 
uniform dark matter starting just from the halos of galax- 
ies, in which case the standard FRW model works. But 
then the fractal distribution of luminous galaxies can ap- 
pear only from special initial perturbations of FRW. n 



2 But, de Vega et al. ( |l99J ) and de Vega & Sanchez (1999) 
showed that Newtonian self-gravitating N-body systems have 
a quasi-equilibrium fractal state with a dimension of ~ 2. 



The second solution is a very low value for the global 
average density (Baryshev et al. 1998; Humphreys et al. 
1998b). However, if the upper cut-off scale of the fractal 
structure is large, the low density contradicts the esti- 
mated density of the baryonic luminous and dark matter. 

4-2. On the applicability of the LTB model to fractals 

The LTB model has proved useful for understanding the 
kinematics of galaxies around individual mass concentra- 
tions. For example, Teerikorpi et al. (1992) could put in 
evidence the expected behaviour in the Virgo superclus- 
ter: 1) Hubble law at large distances, 2) retardation at 
smaller distances, 3) zero-velocity surface, and 4) collaps- 
ing galaxies at still smaller distances. 

Bonner (1972) was the first to apply the LTB model 
to the hierarchical cosmology. He used de Vaucouleurs 's 
density law p ~ eT 7 with 7 = 1.7. Ribeiro (1992a, 1992b, 
1993) has developed a numerical approach to solving LTB 
equation for fractal galaxy distribution. Humpreys et al. 
(1998b) gave a relation between number counts and red- 
shifts for LTB models with large scale FRW behaviour. 

However, the application of LTB models to a frac- 
tal distribution leads to a conceptual problem, because 
the original LTB formulation contained a central point 
of the universe, around which the density distribution is 
isotropic. In a fractal distribution (Mandelbrot 1982] ) there 
is no unique centre, but every object of the structure may 
be treated as a local centre which accommodates the LTB 
centre. Every structure point is surrounded by a spheri- 
cally symmetric (in average) matter distribution. 

In this sense, the application of the LTB model to frac- 
tals means that there is an infinity of LTB exemplars with 
centres on every structure point. Their initial conditions 
are slightly different, because for any fixed scale the aver- 
age density is approximately constant. For different scales 
the density is a power law. This excludes geocentrism and 
makes possible the use of LTB models as an exact general 
relativistic cosmological model where expansion of space 
becomes scale dependent. ^ 

4-3. Simultaneous bang time 

We showed in Sect. 3 that a unique bang time and con- 
stant density imply open FRW models. Here, with non- 
linear LTB models, we study density perturbations with 
arbitrary amplitude on the FRW background and show 
how the initial fractal density changes the models. 



Humphreys et al. (19981) used a density profile which 



needed junction conditions for densities corresponding 



3 A similar property exists for Friedmann models where the 
rate of space expansion within distance d from a fixed galaxy is 
determined by the total mass of the sphere around this galaxy. 
For LTB fractal models the space expansion at distance d from 
a fixed point of the fractal structure is also determined by the 
average mass of the sphere around this point (Paper III) . 
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Fig. 1. The solution of Eq. (p0[) with unique bang time and 
initial conditions. £,zv, the solution of Eq.(30), depends 
on the initial conditions (M) and (p5[). The upper curve 
is the function log (rFRwiio/^*) 1 ^ 2 ) f° r ^o = 0.001 and 
A = 0.02 (see Table 3), while the lower curve corresponds 
to Vlo = 0.001 and A = 0.002 (Table 4). The upper line 
is maxlog(^I /Cl ) = log(arcsin(l)) = 0.196. The lower line 
is min log(^I /Cl ) = log(|) = —0.176. The galaxy scale e 
is 2 • 10~ 6 . The curve intersects the upper line at £tb- 
It is seen that Slo = 0.001 and A — 0.02 produce the 
intersection at a scale > galaxy scale, which corresponds 
to real (and positive) solution of Eq.(p0|), but Slo = 0.001 
and A = 0.002 do not produce it. In the last case the 
solution of Eq. (pOj) is imaginary. 



Fig. 2. The coordinate criterium of existence of £zv f° r 
simultaneous bang time tfrw = 0.98, Qq = 0.01 and 
A = 0.02 (see Table 3). The model is defined for £ > izv ■ 
At the boundary £zv velocity £,{£,zv) — 0, see Fig. 4. 

to different scales. We consider the analytical case of a 
smooth density profile. The fractal density on the FRW 
background and simultaneous bang are the initial condi- 
tions of our LTB models: 



*&0) 



1 



.4 



e + Co 



1. 



(34) 



^(£o) = Tf (oo) = TF RW = COHSt. 



Tr-^U^]/^-^ 1 -^) ^ 



Here e = R 



galaxy / ^0 



/h 



10 Kpc/ 5-10 6 Kpc = 2-10" 6 . Above 
the galactic scale, Eq.(|i|) describes the fractal density law 
with D = 2. The density contrast {p ga iaxy/ Pfrw) is <5(£ = 
0) - 10~ 24 g J cm 3 /I0~ 29 g/ cm 3 = 10 5 . So, A ~ 0.2. For 
our calculations we use A = 0.002, 0.02, 0.2, 2, which im- 
ply the amplitude of the density 6(0) = 10 3 , 10 4 , 10 5 , 10 6 , 
and we use O = 0.001, 0.01, 0.1, 0.99, which imply 
tfrw = 0.997, 0.98, 0.898, 0.688. 

The properties of the LTB models with initial con- 
ditions Eq. (p4]) , Eq.(p5|) were studied in Sect. 3. Now we 
apply the results of Sect. 3.2. to the initial conditions with 
given parameters. For the chosen values of A and Oo LTB 
models have a closed "core" and open "shell". But only 
for = 0.001 and A = 0.002 does Eq.@ have a complex 
solution (Fig.l). This means that only these parameters 
produce the TB model with fractal density and simulta- 
neous bang time with ^zv > (see Table 4). For all other 
cases of the adopted parameter values, Eq.(|30|) has a real 
solution and £zv > (see Tables 1-4). 

Figs. 2 and 3 illustrate the coordinate and mass criteria 
for the existence of £zv > for A = 0.02 and O = 0.01. 

Our solution depends on Slo : M* ~ ^o & n d tfrw — 
tfrw(^o)- Tables 1-4 show characteristic values of 



Fig. 3. The mass criterium of existence of £zv > for 
simultaneous bang time tfrw — 0.98, Q.q = 0.01 and 
A = 0.02 (see Table 3). The upper of two parallel lines 
corresponds to /itb and the lower line corresponds to \ifi , 
see Eq.(|l|) and Eq.@. The model is defined for Co > 
t;zv- At the surface £zv velocity £(£zv) = 0, see Fig. 4. 

£zv and £fi for different A and flo. Here ^o = c /Hq = 
5000 M pc, H = mkm- l s- x Mpc, 

Finally, we calculate the velocity £(£), produced by the 
initial conditions Eq.(B4]) and Eq.(p3). For both domains 
of the model, closed and open, 

(36) 



Zo = VW/So)B, 

where B is the solution of the equation (27). Fig. 4 shows 
the resulting non-linear velocity-distance relations down 
to the zero- velocity radius, which in a natural manner 
marks the distance close to which the deflection from the 
Hubble law becomes very significant. 

Tables 1-4 and Fig. 4 confirm the previous conclusion 
by Baryshev et al. (1998) that the observed linear Hubble 
law is compatible with such a fractal density only if the 
FRW density parameter Oo is small. For instance, if Oo = 
0.99, then the zero-velocity radius Rzv ranges 3.4 - 3400 
Mpc for the range of the density contrast A = 0.002 — 2, 
and is 344 Mpc for the "preferred" value of A = 0.2. With 
a very small value of O , 0.001, Rzv appears around 0.6 
Mpc, which is an intergroup scale, while a good linear 
Hubble flow is reached around 6 Mpc. 

4-4- Nonsimultaneous bang time 

Here we study still another solution of the HdeV para- 
dox within the LTB model: a non-unique bang time may 
produce the linear Hubble law within a fractal structure 
having an upper cut-off. 

We take the initial conditions from the observations in 
the domain from 1 Mpc to 100 Mpc, i.e. the fractal density 
and the Hubble law: 

<5(£,0) = -^- + l, i = t (37) 

e + £o 

The bang time in this situation is calculated by two for- 
mulas, corresponding to the closed and open domains: 
1 



Tf(&) = 



(1-B) 3 / 2 



1 



-VsVT - H 



T ^ 0) ~ (B 1)3/2 



3in y/l-B 



y/By/B- 1- 



arcsm 



hy/B- 1 



(38) 



(39) 



where from Eq.@ and Eq.@, B = £§/ju*(£o)- For A = 
0.02 and O = 0.01 the bang time T£(£o) is shown in Fig. 5. 
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Fig. 4. Velocity produced by initial conditions Eq.(M) 
and Eq.(p3) (simultaneous (FRW) bang time and given 
density profile) for il = 0.001, 0.01, 0.1 (from left to right) 
and A — 0.02. Different £Iq produce different £,zv, in which 
£,(£,zv) = 0. Recall that log£ = — 1 corresponds to 500 
Mpc and the size of the Local Group (1 Mpc) is encoun- 
tered at log£o = —3.7. 
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Table 1. A = 2, 6(0) = 10 6 



Fig. 5. The non-simultaneous bang time T£(£o) produced 
by initial conditions @ with Q Q = 0.01, A = 0.02. At 
infinity the time of collapse t^(^ ) — ► tfrw = 0.98, indi- 
cated by the horizontal line. 



5. Discussion and conclusions 

We have studied the LTB model solution for two pairs of 
initial conditions related to observations: (bang time func- 
tion, density profile) and (density profile, velocity func- 
tion) which both permit one to treat the LTB problem as 
the Cauchy problem. The second pair may be obtained 
from observations at the present epoch, while from the 
first pair, with bang time function assumed, one may pre- 
dict the present-day velocity function. A discus sion of this 
mathematical side is given by Gromov ( 1999 ). We show 
that for most parameter values of the model, there is the 
zero- velocity surface. 

Our study is much guided by Bondi's (1947) idea when 
he said "The assumption of spherical symmetry supplies 
us with a model which lies between the completely homo- 
geneous models of cosmology and the actual universe with 
its irregularities." 

As the "actual universe with its irregularities" has 
turned out to be fractal, at least in its luminous matter 
and in scales up to 100 Mpc or more, we follow some ear- 
lier works in representing fractality with spherically sym- 
metrical systems of different scales. The interesting con- 
ceptual difficulties will be treated in Paper III. Here we 
have assumed that the fractal representation is adequate 



and complement our previous work (Baryshev et al. 1998) 
on the Hubble law within fractals. The qualitative conclu- 
sions of that paper, based on the linear regime, are con- 
firmed by our exact LTB solutions. In particular, a fractal 
distribution of matter with D = 2, smoothly going over 
to the FRW background, generally results in a large devi- 
ation from the linear Hubble law, when there is a unique 
bang time. Only if Slo ~ 0.001, does a reasonable density 
contrast A = 0.2 produce an acceptable Hubble law. 

As the above mentioned two pairs of initial conditions 
are interconnected, one may start from the linear Hubble 
law as the velocity function, and derive the required bang 
time function. This we have done, and conclude that in 
this manner in the frame of the LTB models it is thus 
possible to have a linear velocity - distance relation when 
matter distribution is fractal. Physically, this would mean 



«0 


£,ZV 


Rzv(Mpc) 


*/i 


R f i(Mpc) 


0.001 


1.2 10" 4 


0.61 


6.7 10" 4 


3.3 


0.01 


0.001 


6.1 


0.007 


33 


0.1 


0.01 


53 


0.06 


332 


0.99 


0.07 


344 


45 


210 5 



Table 2. A = 0.2, 5(0) = 10 5 
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Table 3. A = 0.02, 6(0) = 10 4 . 
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Table 4. A = 0.002, 6(0) = 10 3 



that different spherical shells are created at different mo- 
ments given by the bang time function. 

Thus the list of possible solutions of the Hubble-de 



Vaucouleurs paradox (Baryshev et al. 1998) in the frame 
of the LTB description of fractality now includes 1) a very 
low cosmic density, or 2) a dominating smooth dark mat- 
ter, or 3) non-simultaneus bang time. 
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